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Abstract 

We compute B — > Krf^ branching ratio using perturbative QCD ap- 
proach. We show that a triangular relation among amplitudes for B° — > 
K°-k°, B° — > K°r], B° — > K°r]' receives large corrections from SU(3) breaking 
effects. If experimental value will come closer to the lower limit of the present 
BELLE data there will be a possibility to understand the large branching ra- 
tio of B° — > K°r]'. Otherwise, we perhaps need to modify our understanding 
of rj' meson, for example, inclusion of a possible admixture of gluonium state. 
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Three years has passed since CLEO announced an unexpectedly large 
branching ratio for B — > Krj decays H: 

Br(B° -> KS) = (89±i| ± 9) x KT 6 (1) 

BELLE also reported their results in BCP4 conference, 0: 

Br(B° -> i?V) = x 10 ~ 6 ( 2 ) 

Various theoretical suggestions have been made to understand the large 
branching ratio. While new physics contributions were discussed [[| |]] we feel 
that better understanding of the standard model calculation of the branching 
ratio is necessary. In Ref. ||, it was shown that there is a possible choice 
of theoretical parameters involving form factors, CKM parameters, nonfac- 
torizable contribution and decay constants of 77 — rj' system which gives a 
branching ratio consistent with the experimental data. A SU(3) relation 
which is independent of most of the above mentioned uncertainties has been 
derived § % |: 

- 3y/2 A(B° -> K°tt°) + aV3A(B° -> K°rj) = VQA(B -> K°r,'). (3) 

Using the CLEO measurement 0, Br(B° -> A°tt ) = (14.6±|;?±|;|) x KT 6 , 
and theoretical expectation that |A(i? — > i^°?7)| is small compared to the 
other two amplitudes, the observed value for Br(B° — > K°r]') in Eq. (|1]) 
seems to be too large. This relation also excludes explanations which make 
Br(B° — > A' ^ ) increase simultaneously with Br(B° — > K rj'), for instance, 
invoking large Wilson coefficients with new physics effects, or increasing the 
input parameters like form factors, the CKM parameters, etc. 

In this letter, we perform calculation of the branching ratio by using per- 
turbative QCD (pQCD) approach and examine the B —>■ Krj' problem. The 
SU(3) breaking effect is included through the decay constants and the wave 
function and as a result, Eq. (|j) is modified. We also give a theoretical esti- 
mate of color suppressed penguin contributions which is one of the candidate 
mechanism to enhance Br(B° — > K°r)('>), but not Br(B° — > A°7r°). 

In the 80's, rj gluonic admixture was examined in Ref. ||10|| . Recently, 
there have been some progress in understanding the 77 — 77' system. We use 
the 77 — 77' mixing angle and the definition of the decay constant in 77 — 77' 
system which include recent improvements. A simple description of 77 — 77' 
states was introduced in the literature, 
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where X„{,) , Y-w and Z v i parameters represent the ratios of uu + dd, ss and 
gluonium component of rf'\ respectively. This work was updated by one of 
the authors [II]. In this work, the gluonium content of rj' is reanalyzed using 
all available radiative light meson decays and a result 



Xtf + Y„> 



< 0.26 



(6) 



is obtained, which indicates that 26% of gluonic admixture in rf is still pos- 
sible. 

Ignoring the small tree contribution, we can write the amplitudes of B° — > 
K°tt° and B° -> K°rj^ decays as: 



A(B°^K%) = -l/V2P d 

A(B° — > K°r}) = X v /V2P d + Y n P s + P 

A(B°^K°rj') = X v ,/V2P d + Y v ,P s + P', 



(7) 
(8) 
(9) 



where Pd(s) includes color allowed bsdd(ss) penguin and annihilation pen- 
guin contributions and P^ is SU(3) singlet contribution. We depict the 
corresponding diagrams in Fig.|l]. 
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Figure 1: Diagrams for bsdd penguin, P d , bsss penguin, P s and SU(3) 
singlet penguin, P' contributions 

Here we introduce a parameter r which represents the SU(3) breaking 
effect: 

P s = rP d (10) 
and a parameter which represents the ratio between and P s : 



P = sR 



P' = s'P s . 
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Using these parameters, Eqs. (0-0) lead to £77(3) relation in general form 

A(B° -> K°r)) = -{X v + rV2Y v + rV2s)A(B° -> K°ir°) (12) 
A^^i^V) = ~{ X n> + rV2Y v , + rV2s')A(B° -> ATV). (13) 

Now, let us examine what it takes to obtain Eq. (H). The following 
assumptions must be applied: 

• rj' does not have the gluonic content (Z^ = 0) so that X v r v n and Y v ( v ^ 
are related to the pseudoscalar mixing angles as 




cos a p — sin a p 
sin a p cos a p 




(14) 




For convenience, we also display r\ and rf states in terms of singlet state 
r/i and octet state rj s as 

(15) 

where r)$ = (uu + dd — 2ss) / \/o and 771 = (uu + dd + ss) / a/3. a p can be 
written in terms of the pseudoscalar mixing angle 9 P as a p = 9 P — 61 + 1 
with the ideal mixing 9j = arctan(l/y / 2). The allowed range of values 
for 9 P is —20° to —10°. To obtain Eq. (§), the angle a p is fixed at 
cos a n = V2/y/3 and sina p = l/y/3, which corresponds to 9 p ~ —19.4°. 
We should remind the reader that we are taking a value of 9 p which is 
at the edge of the allowed region. While we stretched the error bar to 
allow the region 9 P ~ —20°, the experimental data for to — > ryy decay 
and 77 — > 77 decay disfavor this region. The best fit range is: —17° to 
-10° in 0. 

• SU(3) symmetry is exact so that Pd = P s , i.e. r = 1. 

• The ratio of the SU(3) flavor singlet contribution to B — > Kr] and 
B — > Kr] 1 is written as s'/s = —cos9 p / sin 6^, which is extracted from 
the ratio of the SU (3) singlet component of r\ and 7/ states in Eq. (|15|). 
It can be easily seen below that this is valid only if we set fx = fn and 
ignore the electric penguin correction factor (£ = 1). 

We would like to point out that whether the experimental value Eq.([I|) is in- 
consistent with Eq. (HD depends crucially on the assumptions above. Looking 
at the SU(3) relation in general form in Eq. (|T3"D, we see that the ampli- 
tude of B° — > K°7]' can be enhanced by large r and s'. And in fact, re- 
laxing above assumptions, we can easily have Br(B° — > K°rf) consistent 
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with data. For example, if we choose r = 1.1 and 6 P = —10°, keep the 
relation s'/s = —cos6 p / sin6> p and take Br(B° — > K 7y ) = 15 x 1CT 6 and 
Br(B° -> K°T]) = 0, Eqs. (TTJ) and (|TJ) give 

Br(B° -> fCV) = 84 x HT 6 . (16) 

We insist that SU(3) breaking effects for r, s and s' must be studied before 
we conclude that experimental data is too large - and that we need new 
physics to explain the observations. The large value in Eq. ([16]) is due to the 
fact that constraints: Br(B° — > K°i]) = and s'/s = —cos8 p / sin# p leads 
to large s'. 

In our pQCD approach, we evaluate these SU(3) breaking parameters as 
well as the branching ratios Br(B° — > K°r]') and Br(B° — > K°7]). Now let us 
explain our calculation for B° — > K°r)^ decay amplitude. (Full calculation 
will be presented in elsewhere.) The pQCD approach is developed to give 
more precise theoretical prediction beyond vacuum saturation approximation 
13 1 . In this approach, the amplitude for B° — > K°r]^ is given as (see 



also Fig. g): 

M(B^K^) = ^{VJ K ^Fi + V t f y Y v( , )F » 

+V t {f x ^ + f x ^ + f y Y v{l) )Fr 

-VJ^F 1 / + Vth^Fi + V t f B Y v(l) F a B }, (17) 

where £ is a correction factor which represents the difference between electric 
penguin contributions of bsuu and bsdd(bsss) penguin diagrams and V q = 
V* s Vqh (q = t,u). The parameter aj(/i) in Fig. || is defined as a2(/i) = 
Ci(fi) +C 2 (fi)/N, ai(fi) = Ci(fi) + C i+ i(n) /N for i = 3,5,7,9 and cti(/i) = 
Ci(fi) + Ci-i(n)/N for i = 4,6,8, 10, where N is number of the color and 
Ci(fi) is Wilson coefficient (we use the same definition of Wilson coefficient 
as the one in |12|). It is worthwhile pointing out that scale fi in these 
coefficients is related to the loop integration variable for diagrams shown 
in Fig. 2. Thus the usual scale dependence problem associated with the 
factorization assumption is absent in the pQCD approach. 

In our notation, F\ (i — I ~ III(IV)) represents penguin (tree) diagram 
and F l a [i = A, B) represents annihilation penguin diagram. The correspon- 
dence between Eq. (|j) and (|17|) is as follows: 

Pa = ^-VMkFI + f B F?) (18) 
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Figure 2: Contributing diagrams. 
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Note that we also have additional contribution from nonfactorizable diagrams 
(see, Fig. |^), which is not calculated in vacuum saturation approximation 
but can be calculable in pQCD approach. However, we found that nonfac- 
torizable contributions to branching ratios for B° —>■ K°1]' and B° — > K°r] 
are less than 10%. 




Figure 3: Nonfactorizable contribution 

As is mentioned in the introduction, we use new definitions of decay 
constants in the 77 — 77' system - the decay constants at the non-anomaly limit 
BHU: 



i/«P„ = <0\(ur l5 u + dr l5 d)/V2\^J^> 



(21) 



ijyPp. = < 0\sYlbs\ss > . (22) 



For the value of f x and f y , we use the ones that are given in |L4[ where isospin 
symmetry is assumed for f x and SU(3) breaking effect is included for f y : 



/. = /*, fy = pP K -fl (23) 

These values are translated to the values in the two mixing angle method, 
which is often used in vacuum saturation approach as: 

/ 8 = 169 MeV, f x = 151 MeV, (24) 

# 8 = -28.9°(-18.9°), 0i = -10.1°(-0.1°), (25) 

where the pseudoscalar mixing angle 9 P is taken as —20° (—10°). The wave 
function for dd components in r/^ and K meson are given as: 



%JP, x, C) = frfjx) + ™ V o d *rfJx) + Cn#*(M ~ v ■ n)r v 'Jx) (26) 
* K (P, x, C) = f (j>£(x) + m^^(x) + Cm K (# - v ■ n)<f>%(x), (27) 



where P and x are the momentum and the momentum fraction of r] dd -(K), 
respectively. We assumed here that the wave function of rj dd is same as 
the 7r wave function. The parameter ( is either +1 or —1 depending on 
the assignment of the momentum fraction x. 4>^ d3 (K)' ^ -(K) anc ^ -(K) 
represent the axial vector, pseudoscalar and tensor components of the wave 
function, respectively, for which we utilize the result from the Light-Cone 
sum rule |17J including twist-3 contribution: 

O) = ^/,x(l-x)[l + a^(5(l-2a;) 2 -l)] + 

af*— (21(1 - 2xf - 14(1 - 2x) 2 + 1)] 
8 

= 2^ / ' [1 + (3 ^"§^ (3(1 " 2a:)a " 1) + 
97 81 1 

(- 3 ^3 - Yo P 2 r ldd - ^^(35(1 - 2xf - 30(1 - 2xf + 3)] 



<p a '(x) = -^=fjl-2x) 



[\ + (5^3 - |w - Y/ ndd - - \p 2 Vu af° O(10* 2 - 10x + 1)] 
<f£(x) = -^=f K x(l-x)[l + 3a[ < (l-2x) + ^(5(l-2x) 2 -l)], 
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where 



af* = 0.44, af 3 = 0.25, of = 0.20, of = 0.25, 
^w = ^&. % = 0.015, ^ = -3. 



i) 



We assume that the wave function of uu is same as the wave function of dd. 
For the wave function of the ss components, we also use the same form as dd 
but with m s s and f y instead of and f x , respectively. The pseudoscalar 
and tensor components of the K wave function are obtained by exchanging 
parameters of the pseudoscalar and tensor components of dd wave function, 
respectively as follows: 

fx <-> fx, P Vdd <-> Pk = 5' a 2(i a iV (28) 

The parameters m % Q (i = rj^uu) , r} a s, K) are defined as: 



m ' = m = ■ -, = — - — r — , m 



(m u + m d y "" u " (2m.) ' u ~ m d(u) + m/ 

(29) 

Because there are large ambiguities in quark masses, parameters defined in 
Eq. (p9|) introduce considerable theoretical uncertainties. In our analysis, we 
use constraints for these parameters from analysis of other decay channels. 



In Ref. [0 branching ratios for B — > ixix are analysed in pQCD approach. 
The allowed region for is given as l.lGeV < < 1.9GeV. Ref. |T2| , 
which studies B — > Kir in pQCD, gives the best fit value of m s = 140 MeV. 

We saw that the SU (3) breaking effects were included through the decay 
constants and the wave functions in Eq. ( p3| ) and Eq. (|29|), respectively. In 
exact SU(3) symmetry limit, 

fn = fx = fy = fx, (30) 

mf^ = mf = m", (31) 

Eq. (|[)is recovered. 

Now we show our numerical results. The parameters which are used in 
our calculation are as follows: 

G F = 1.16639 x 10~ 5 GeV" 2 , A ^ = 250 MeV,a s (M z ) = 0.117, a em = 1/129, 
r B o = 1.56ps, A = 0.2196, A = 0.819, R b = Vp 2 + V 2 = 0-38, 3 = 90°, 
M w = 80.2 GeV, M B = 5.28 GeV, m t = 170 GeV, m b = 4.8 GeV, 
f B = 190 MeV, f K = 160 MeV, f w = 130 MeV 



7 





-1.57 (-1.93) 


/^'^ 


-4.03 (-3.50) 


(UX v ,/V2t+f x X v ,/V2 


2.16 (2.15) 




i0.0334 (z0.0411) 




0.133 (0.163) 




0.734 (0.902) 


f B Y v rReFB 


0.494 (0.428) 


f^lmF* 


2.00 (1.74) 





-2.26 (-1.95) 


f Y F H 

Jy V e 


2.81 (3.50) 


4-f Y \F nI 


0.120 (-0.266) 


Vu Xr, pIV 

Vt V2 r e 


i0.0481 (z0.0417) 


f B ^ReF a A 


0.191 (0.165) 




1.06 (0.915) 


f B Y v ReFB 


-0.343 (-0.428) 


f B Y,lmF* 


-1.39 (-1.74) 



Table 1: 



The numerical result with the best fit parameter set, 



{mo* s, ,m 8 ,w B } = {1.4GeV,140MeV,0.4GeV} for the each term in Eq. 
(P~7|) . The pseudoscalar mixing angle is taken as 9 P = —20° (—10°). The 
values are factored out by 10 4 Vj. 



The B meson wave function is given as follows: 

<Mx) = N B x\l - x fex V \-\{^f - 
N B = 91.7835 GeV, 

where uj b is a free parameter. At first, we show our result with the best 
fit parameter set which is obtained from analysis of B — ► tttt and B — > Kn 
processes in pQCD approach [O, III 



{m v dd(uU \m s ,u B } = {1.4GeV, 140MeV, 0.4GeV} (32) 

Our result of the branching ratio for the above parameter set is 

Br(B° -> K°T]') = 18(16) x 10" 6 (33) 
Br(B° -> K°rj) = 0.44(1.9) x 10" 6 , (34) 

for the pseudoscalar mixing angle 9 p = —20° (—10°). The computed am- 
plitude for the each diagram in Fig. |2] is given in Table [1|. With the 
same set of parameters, we obtain theoretical prediction for B° — > K°7r° 
as Br(B° K°tt°) = 9.5 x 10" 6 . 

Since we observe large imaginary part for annihilation diagrams we need 
to rewrite the SU(3) relation in Eqs. ( |T2"D and ( |T3"D more precisely. Again 
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ignoring the small amount of tree contribution, we obtain: 



ReA{B° - 


^K° V ) = 


-{X v + rV2Y v + rV2s)ReA{B° -> K°ir°) 


ReA(B° - 


+ K°r)') = 


-{X rf + rV2Y rf + rV2s')ReA(B° -> K°ir° 


lmA(B° - 


^K°rj) = 


-{X v + r'V2Y v )lmA{B° -> K°ir°) 


lmA(B° - 


+ K°rf) = 


-(X v , + r 1 V2Y v ,)lmA(B° -> ifV) 



(35) 
(36) 
(37) 
(38) 



where 



fFcFi+fsReFA- 
fyFll+f B ReFB 



J _ JmF n B 

fyFl'+fBReFB 



(39) 
(40) 



Using the values listed in Table [l], the SU(3) breaking effect r^'\ which is 
assumed as r = 1 in Eq.(^), is calculated as r = 1.2 and r' = 1.3 and the value 
of s(> which is proportional to the SU(3) singlet contribution is calculated 



as s 



-0.50(-0.50) and s = -0.028(0.059) for 9 V 



-20° -10 c 



The 



electric penguin correction factor is obtained as £ = 0.543. To employ the 
new definition of the decay constant in 77 — rj system modifies the relation 
between s and s' as 



s / s 



V2 



(41) 



which is assumed to be s' /s — — cos Q p / sin 9 P in Eq. (|3]). s' always has minus 
sign. This is due to the sign difference between F£ and F£ , which can be 
traced back to the relative size of the Wilson coefficients. This effect is also 
seen in calculations using vacuum saturation approximation [|l^]. This fact 
implies that the SU(3) singlet penguin contribution to B° — > K°r]' tends to 
decrease the branching ratio. 

Before we show our numerical results for the different parameter sets of 
the wave functions, it might be convenient to summarize the trend of the size 
of the amplitudes for the variation of the parameters, diagram by diagram in 



Fig. |. 



ml • 11 j r ldd(uu) 

1 he variable parameters are m 



and rriQ which depend on 



quark masses and lob which parameterizes momentum distribution of b quark 
in a B meson. The amplitudes Fj, ReF A and lmF A depend on parameters 
d m^, Fg 1 , RgFq and ImF^ depend on parameters m[j ss and 
and F£ depends on uiq . All the amplitudes are increased when ujb is 



Vdd(uu) 

m 



K 



m 

decreased. Note that the branching ratio of 5 C 
of F/, ReF a A and ImF a A . 



K 7r are written in terms 
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We first discuss the input parameter dependence of s' and s. When we 
decrease m s , which increases mj} ss and ttIq , the amplitudes F^ 1 , ReF^ and 
Fl u get enhanced simultaneously. When we decrease ub, F^ 1 and F^ 11 
increase while Rei 7 ^ remains unchanged. Therefore, as we can see from Eq. 
d40|) , s' and s are quite insensitive to the input parameters mo", and ujb- 

Now let us see how we can obtain large Br(B° — > K°i]'). To enhance 
Br(B° — > K°T)'), we need large r. The value for r is increased for smaller 
mQ dd(uu) , smaller m s or smaller ujb- First, let us try lub = 0.3GeV which is 
the lower limit considering other two body B decays. For example, we obtain 

r = 1.2, r' = 1.3, s = -0.032(0.054), s' = -0.49(-0.49) 
Br(B° -> K°n°) = 21 x 10~ 6 , Br(B° -»■ J? //') = 34(31) x lO" 6 , 

for 9p = — 20°(— 10°). In order to further increase r and thus increase B° — > 
K°T}', we need to reduce m s , however, it ends up also increasing B° — > 
K°7r° branching ratio. This is unacceptable. Therefore, we put back ujb as 
0.4GeV and try to enhance r only by changing ml dd< - uu) and m s . In fact, the 
dependence of r on mQ dd(uu) is so weak that even if we take smaller value for 
m^'"™' the branching ratio for B° — > K°r]' does not become large enough. 
So, the only possibility is to take a smaller value for m s . The results evaluated 
with {m s ,uj} = {100MeV,0.4GeV} are as follows: 

For m^" 81 = 1.9 GeV 

r = 1.4, r' = 1.3, s = -0.025(0.045), s' = -0.40(-0.40), 

For rn dd(uu) = 1.4 GeV 



r = 1.7, r' = 1.5, s = -0.025(0.045), s' = -0.40(-0.40), 



For ml muu) = 1.1 GeV 

r = 2.0, r' = 1.6, s = -0.025(0.045), s' = -0.40(-0.40), 

Obtained branching ratios for B° — > K rj' and B° — > K°ir are given in 
Tablep]. Considering that BELLE reported a smaller branching ratio for 
B° — > K rj' and also that experimental value of B° — > K°7r° has still large 
error and can be large, a situation such that mQ dd{uu) = 1.9 GeV can not be 
excluded. 

For a comparison to other approaches, we give the obtained values of the 
form factors for different parameter; 

F^{0) = (0.36, 0.30, 0.26) : rn dSiu9) = (1.9, 1.4, l.l)GeV 
F^-* K (0) = (0.47, 0.35) : m s = (100, U0)MeV 

where u>b = 0.4GeV. 
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Br(B° -> A°7r°) 


Br(B° -> A<V) 


Br(B° -> A ?]) 




21 x 10" 6 


50(45) x 10" 6 


2.4(7.5) x 10" 6 


m ^d(uu) = lA GeV 


15 x 10~ 6 


45(39) x 10" 6 


4.6(11) x 10" 6 




11 x 1(T 6 


41(35) x 10- 6 


7.0(13) x 10- 6 



Table 2: The numerical results for the branching ratios of B° — > A°7r°, 
B° — ► A ?]' and 5° — > A ?] with parameters tuned to increase 5° — > K°rj' . 
{m s ,u} = {lOOMeV, 0.4GeV} and different m n dd{uu) . These parameter sets 
are allowed by the pQCD analysis for B — > mr and B — + Kti processes. The 
pseudoscalar mixing angle is taken as 9 P = —20° (—10°). 



In fact, there is another interesting aspect. Looking at Eq. dl2|) , the first 
and the second term which are the dominant contributions have an opposite 
sign so that the relative variation of the branching ratio for B° —>■ K°r) is 
much larger than that for B° —>■ K°i]' when r varies. The branching ratio 
for B° — > K°r], which is considered to be negligible in Eq. @, is enhanced 
greatly depending on the parameters. We summarize our numerical results 
of the branching ratios for B° — > A°7r°, B° — > A°r/ and _B° — > A°?7 in Table 
0. We can see a large dependence of the branching ratio for B° — > A°?7 on 

As is mentioned before, # p = —20° is the smallest limit of the allowed 
region. Hence, our results for 8 P = —10° indicate that B° — > K°rj process 
may be observed soon. 

In conclusion, we examined the large branching ratio of B — > Krj' process 
using SU(3) relation in general form, Eqs. (|T2"| ) and (p!3|). If there is a large 
SU(3) breaking effect, which means that r is much larger than 1, or there is 
a large SU(3) singlet penguin contribution, which means s' and s are very 
large, Eqs. and (|13D imply that we would have large branching ratio 

for B° — > A°r]'. We computed r, s' and s as well as branching ratios for 
B° — > A°7r°, 5° — > K°i]' and 5° — > A°r] processes in pQCD approach, s' 
is found to contribute destructively to the other dominant contributions to 
B° — > A°r]' process. Our numerical result in Table |2| indicates that in a case 
that the experimental data for B° — > A ?]' will come close to the lower limit 
of BELLE data, B° — ► A //' problem can be understood in the standard 
model. However, in this case, the correlation of the experimental data for 
B° — > A //' to the experimental data for B° — > A°7r° and B° — > A ?] has 
to be examined carefully. In particular, considering that the relatively large 
value of 6p which is close to —10° is favored by recent experiments, we have 
to keep in mind that Br(B° K rj) may not be so small. If Br(B° — > A°r/) 
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remains high at its present value or the combination of the experimental data 
for B° — > K°n°, B° —>■ K°1]' and B° — > K°rj deviate from our result in Table 
0, it may imply that we need modify our understanding of rf . 

Finally, we would like to make a comment on two suggested mechanisms 
to explain the large branching ratio for B — > Krf, 

(1) Intrinsic charm contribution: 

The Cabbibo allowed b — > ucc process can contribute to B — > Krj' if 
there is intrinsic cc content in 77' (see Fig. [|(a)) (1^, 20]. 



The amount of the cc content in rf which is parameterized by decay constant 
fv' as: 

iftv* = (0\c^ 5 c\rj'(p)) (45) 

Decay constant f$ is obtained from the radiative J / ip decay and two photon 
process of v( . It was found that the numerical result for the intrinsic charm 
contribution to B — > Krj' decay was very small [13, [15 . 



(2) The SU(3) singlet contributions: 

B — ► Kv( is produced by fusion of gluons, one gluon from b 
process and another one from spectator, (see Fig. 3(b)) pT] , 22. 23 . 





Figure 4: (a) Intrinsic charm contribution, (b) SU(3) singlet contribution 



According to the paper [IT], there is still a possibility that 7/ includes at 
most 26 % of pure gluonic state, gluonium. The contribution of the diagram 
in which two gluons in Fig. |](b) are directly attached to gluonium in rj' 
instead of attached to triangle quark loop may be important for B — > Krj' 
decay. 
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Note added in proof 

After this paper was submitted for publication, BELLE and BABAR an- 
nounced new data: 





Br(B° -> K"i) 


Br{B+ -> HT+jj') 


BELLE |25 


(551$ ± 8) x 10- s 


(791^ ± 9) x 10~ B 


BABAR JZB| 


(42-11 ± 4) x 10~ B 


(70±8±5) x 10~ e 
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